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On  Lacunary  Multiresolution  Methods 
of  Approximation  in  Hilbert  Spaces 


Lubomir  T.  Dechevski  and  Wolfgang  L.  Wendland 


Abstract.  We  study  lacunary  multiresolution  methods  from  the  point 
of  view  of  their  analogy  to  the  use  of  near-degenerate  elements  in  finite 
and  boundary  element  methods.  The  main  results  are  characterization  of 
the  best  N-term  approximation  of  solutions  of  nonlinear  operator  equations 
and  best  N-term  approximation  by  near-degenerate  normal  approximating 
families  in  Hilbert  spaces. 


§1.  Introduction 

This  communication  is  part  of  a sequence  of  papers  exploring  the  use  of  near- 
degenerate elements  in  finite-  and  boundary-element  methods  (see  also  [5,6] 
and  their  wavelet-based  analogues,  lacunary  multiresolution  methods.  The 
use  of  near-degenerate  and  lacunary  methods  for  solving  operator  equations 
is  of  considerable  practical  significance  because  in  many  important  problems 
arising  in  industry,  engineering  and  natural  sciences,  the  use  of  such  methods 
leads  to  a dramatic  reduction  of  execution  time  and/or  computer  resources. 
The  theoretical  justification  for  the  use  of  such  methods  is,  however,  very 
challenging:  it  has  been  successfully  carried  out  only  in  a number  of  special 
cases,  by  specific  techniques  which  vary  from  case  to  case.  The  purpose  of 
this  sequence  of  papers  is  to  develop  a general  approach  to  overcoming  the 
challenges  of  the  use  of  lacunary  multiresolution  and  near-degenerate  finite 
and  boundary  element  methods.  Because  of  the  limited  space  available,  we 
shall  consider  only  multiresolution  methods  for  operator  equations  in  Hilbert 
spaces,  with  an  outline  of  the  main  ideas  of  the  proofs,  which  in  the  Hilbert- 
space  case  are  simpler  and  relatively  short.  A much  more  technically  involved 
and  detailed  discussion  of  both  near-degenerate  finite  elements  and  lacunary 
multiresolution  methods  and  the  important  parallel  between  them  will  be 
given  for  more  general  types  of  nonlinear  operators  in  quasi-Banach  spaces  in 
a later  paper. 
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§2.  Approximate  Solutions  of  Nonlinear  Operator  Equations 

In  this  section  we  consider  a general  class  of  nonlinear  operator  equations, 
and  study  the  numerical  solutions  of  these  equations  obtained  by  iterative 
and  projection  methods. 

Let  X,  Y be  real  Hilbert  spaces.  The  class  of  nonlinear  operators  to  be 
considered  is  the  space  LH(X , Y)  of  all  Lipschitz  homeomorphisms  F between 
X and  Y,  that  is,  3 F~ 1 on  Y and  3 C(F,X,Y)  < oo  : ||F(xi)  — F(x2)\\y  < 
C\\xi  — x2\\x,  € X , and  analogously  for  F-1. 

Let  if  be  a Hilbert  space,  such  that  X fi  Y C H and  X fl  Y is  dense  on 
H,  and  let  Y be  the  dual  of  X pivotal  to  H,  i.e.,  the  dual  with  respect  to  the 
duality  functional  defined  by  the  scalar  product  of  H.  We  shall  denote  this 
dual  by  Y = X*  = X*(H).  ' 

Definition  1.  Let  Y — X*(H).  The  (generally  nonlinear)  operator  F : X — > 
X * is  called  Lipschitz,  if 

3C(F,X,H)  < co  : \\F(Xl)  - F(x2)\\x.  < C\\Xl  - x2\\x,  (1) 

VxiVx2  E X,  and  strongly  monotone,  if 

3 c(F,X,H)  > 0 : (F(x i)  - F(x2),xi  - x2)H  > c||xx  - x2|lx>  (2) 

Vxiix2  E X.  The  class  LSM  = LSM(X,H)  consists  of  exactly  those  F : 
X — > X*(H)  which  satisfy  (1,2). 

It  can  be  shown  that  the  constants  C and  c in  (1,2)  are  related  by  c < C. 
It  should  be  noted  that  the  typical  case  here  is  X «— > H <—*  X*  or  I h 
H -f-3  X* , where,  as  usual,  A > B or  B ^ 1 A denotes  continuous  embedding: 
AcFand  ||.||B  <(7||.|U. 

Theorem  1.  (Generalization  of  Theorem  18.5  in  [11]  and  strengthening  of 
Theorem  18.5  in  [15]  for  the  case  of  Lipschitz  operators  in  Hilbert  spaces.)  Let 
X and  H be  Hilbert  spaces  with  the  same  cardinality.  Then,  LSM (X,  H)  C 
LH(X,X*(H)). 

Proof:  (Outline.)  By  duality  arguments,  it  can  be  shown  that  the  cardi- 
nality of  X*(H)  is  equal  to  that  of  X and  H.  Therefore,  since  all  spaces 
are  Hilbertian  with  the  same  cardinality,  there  exist  linear  invertible  oper- 
ators R : H —y  X and  S : H — > X*(H)  which  are  isometric  together 
with  their  inverses.  Hence,  the  equation  F(x)  — y,  x E X,  y E X*(H)  is 
equivalent  to  the  equation  Av  = w,  v,w  E H,  where  Av  = S~lFR.  Now, 
since  F E LSM(X,  H),  it  follows  from  ||5_1||x*—/r  = | |i?| |//— = 1,  that 
A E LSM(H,  H).  Therefore,  by  Theorem  A (see  below),  F is  bijective  from  X 
to  X* . By  a condition  of  the  theorem,  F is  Lipschitz;  it  remains  to  prove  the 
same  for  F -1.  Indeed,  by  the  strong  monotonicity  of  F,  setting  x\  = F~l(y\), 
x2  = F~1(y2),  Vj/iVy2  E X* , we  get 

HF-1^)  - F-1(y2)||x  < \{yi  - 2/2,  F~1(yi)  - F~\y2))H 

< “lll/i  — 2/2 1 |x-  ||F_1(2/i)  - F-1(y2)\\x, 
which  completes  the  proof.  □ 
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Theorem  A.  (See  [11],  Theorem  18.5).  Assume  that  X = X*  = H.  Then, 
LSM(H,H)  C LH(H,H)  and  the  operator  T£<y(x)  = x — e[F(x)  — y],  x £ H, 
is  contractive  in  H for  0 < e < gt,  uniformly  in  y E H,  where  C and  c are 
defined  in  (1,2).  The  best  contraction  factor  is  1 — c2/C2  and  is  achieved  for 
e = c/C2. 

Following  the  idea  of  the  proof  of  Theorem  1,  Theorem  A can  be  modified 
for  the  case  when  X/I.  We  omit  the  details. 

In  the  remaining  part  of  this  section  we  shall  consider  methods  for  ap- 
proximate solution  of  the  equation  F(x ) = y,  x € X,  y &Y,  where  X,  Y are 
Hilbert  spaces. 

Definition  2.  ( See  [12]).  Let  X be  a Hilbert  space.  G C X is  called  an 
existence  set  for  X,  ifVx  e X 3gx  £ G \ \\x  - gx\\x  = min9£G \\x  - g\\x  = 
Eg(x)x  ■ (The  best  approximation  gx  need  not  be  necessarily  unique.)  The 
sequence  {Gat}“=1,  Gn  C X is  called  a normal  approximating  family  in  X,  if 
for  any  N £ IN,  Gn  is  an  existence  set,  with  Gn  C G n+i  and  Gn  — Gn- i C 

Gin- 

Obviously,  an  existence  set  in  X is  closed  in  X (typical  example:  any 
finite-dimensional  subspace  of  A). 

Definition  3.  Let  X be  a Hilbert  space.  The  sequence  {Gn)n= i : Qn  C X 
is  said  to  have  the  strong  approximation  property  (SAP,  for  short)  if  Uw=i  Qn 
is  dense  on  X in  the  inner-product  topology  of  X. 

Let  us  consider  now  the  Galerkin-Petrov  projection  methods.  Let  Pn  : 
X — > X,  Qn  : X — > X be  projectors  with  dim  Pn{%)  = dim  Qn{X)  = N, 
and  PnPn+i  = Pn+iPn  = Pn,  QnQn+ i = Qn+iQn  = Qn- 

Example  1.  ( Galerkin-Petrov  method  for  monotone  operators .)  For  a Hilbert 
space  X,  let  Y = X*(H).  The  equation  F(x)  = y,  x € X,  y e X* , is  replaced 
by  Q*nF(Pn x)  = Q*Ny,  where  Q*N  : X*  — > A*, dim  Q*N(X*)  = N,  is  the 
Banach  adjoint  of  Qn-  The  N x N nonlinear  system  is  determined  by 

(Q*N^(PNx),QNh)H  = {QnVtQn^H-  (3) 

By  Lemma  23.1  in  [15],  it  follows  that  if  F € LMS( A,  H),  where  A and 
H are  separable,  then  (3)  has  a unique  solution  for  N large  enough. 

In  the  case  A = A*  = H,  if  N is  large  enough,  so  that  Q*nF(PnH ) = 
Q*nH  holds,  then,  by  Theorem  A,  (3)  can  be  computed  by  quickly  converging 
contractive  iterations.  For  small  N , the  condition  Q*nF(PnH)  = Q*NH  may 
fail  even  if  F is  linear  and  Pn  = Qn  (see  [2],  Theorem  10.1.1). 

If  F is  twice  Gateau-differentiable,  then  Newton’s  method  can  be  used 
where  the  inverse  matrix  involved  in  each  iteration  is  usually  sparse.  In  gen- 
eral, this  method  needs  an  appropriate  initial  approximation  xq  to  the  solution 
of  F(x)  = y,  but  if  F is  strongly  monotone  and  potential,  that  is,  if  there  exists 
a real  functional  / : A — > 1R  such  that  F = grad  /,  then,  by  Theorem  5.1  in 
[15] , / is  strictly  convex  and  the  solution  of  (3)  is  equivalent  to  minimizing  the 
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three  times  Gateau-differentiable  functional  /.  Hence,  Newton’s  method  con- 
verges to  the  solution  of  (3)  for  any  xq  £ X,  the  rate  of  convergence  depending 
on  the  constant  c in  (2).  This  technique  is  still  numerically  efficient  if  F is  only 
Lipschitz,  and  Newton’s  method  or  its  various  modifications  be  replaced  by 
the  respective  variants  of  the  more  general  F.  Clarke’s  subdifferential  method. 
In  the  case  of  potential  F,  the  Bubnov-Galerkin  method  ( P 'n  = Qn)  coincides 
with  the  Ritz  method  for  minimization  of  /. 

For  projection  methods  (see  Example  1),  the  strong  approximation  prop- 
erty can  be  written  as  lim/v_00  ||(/x  — FV)x[|x  = 0.  A typical  example 
when  Gn  = PnX  forms  a NAF  having  the  SAP  is  when  Pjv  is  obtained  by 
multiresolution. 

By  Theorem  23.3  in  [15],  if  X is  separable,  then  Gn  — PnX,  as  defined 
in  Example  1,  has  the  SAP;  by  Lemma  23.1  in  [15]  and  in  view  of  F £ 
LMS(X,H),  the  solution'  xn  of  (3)  exists  for  N large  enough  and  ||a;jv  — 
x||x  — > 0,  where  x is  the  solution  of  F(x)  — y,  y £ X* . 

Theorem  B.  (Ce a’s  lemma  for  nonlinear  operators  (see  [13],  Lemma  2.8; 
[11],  Theorems  4.1  and  18.8.)  Under  the  assumptions  of  Example  1,  for  F £ 
LSM(X,  H),  let  x = F~\y)  £ X be  the  solution  of  F(x)  = y£  X*(H),  and 
xn  £ X be  the  solution  of  (3).  Then, 

3 C(F,X,H)  < oo  : EN(F~\y))x  < ||F_1(j/)  - *N||X  < CEN(F~\y))x. 


This  result  shows  that  Galerkin-Petrov  methods  (of  any  type  - finite  el- 
ement or  wavelet)  achieves  the  best  approximation  rates  up  to  a constant 
factor. 


§3.  Best  N- term  Approximation 

For  the  general  paradigm  of  best  A-term  approximation  (BNTAP)  we  refer 
to  [12],  section  3.5,  and  [8]. 

Definition  4.  Let  Xj,  Yj,  j = 0,1  be  Hilbert  spaces,  Xi  c-+  Xq,  Y\  c— > Yq, 
and  let  F £ LH(X0,Y0)  n LH(Xl,Yl).  The  NAF  {Gx}x=i  : gn  C Xu  is 
called  near-degenerate  of  order  (A;  a,/?),  A > 0,  a > 0,  0 > 0,  if  it  satisfies  a 
direct  inequality  of  the  type 

3C  < oo  : EN(F-l(y))Xo<C^F  ^)IIXl , Yy£Y1}  (4) 

where  C = C(N),  with  C x Na;  and  an  inverse  inequality  of  the  type 

3D  < oo  : ||x||Xl  < £WA|M|Xo,  i € GN,  (5) 

where  D = D(N),  with  D x N& . The  partial  case  a = 0 = 0 corresponds  to 
a non-degenerate  (regular)  NAF. 
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^(*o)  :=  {f  eX o : I|/|Uj(xo)  = (||/||9Xo  + £[ ?'E»{f)Xo]')1,q  < c»} 

3=0 

(6) 

and  the  real  interpolation  space 

(y0) :={/eX0:  ll/ll(y„,vi)M  = 

(H/llxo  + jt[2jBK(2-i>f-,Y0,Y1)]‘i)1/q  < oo},  (?)) 

3=0 

where  K(t,  f\Yo,Y\)  is  Peetre’s  /^-functional  (see  [2,12]),  0 < t < oo,  s > 
0,  0 < < 1,  0 < g < oo  (with  the  usual  sup-modification  in  (6,7)  for 

q = oo).  (Recall  that  X\  Xq,  which  explains  the  presence  of  the  saturation 
term  ||/||^o  in  (6,7).) 

Theorem  2.  (Characterization  of  the  best  N-term  approximation  of  solu- 
tions of  nonlinear  operator  equations  by  near-degenerate  NAF  in  Hilbert 
spaces).  Assume  that  the  conditions  of  Definition  4 hold.  Let  0 < q < oo. 
Then, 

(i)  if  0 < a < A and  s : 0 < s < A — a,  then,  3Ci  < oo  : 

HiJ,-1(j/)IU;(x„)  < C1[||F-1(0)||x1  + ||y||(y0lyl)Ti_ii];  (8) 

(ii)  if  (3  > 0 and  0 < s < A + (3,  then  3(72  < oo  : 

IMkvo.yi)^,,  < C,2[||F(0)||y,  + ||F_1(y)m?(x0)].  (9) 

Proof:  (Outline.)  By  a standard  technique,  typical  for  BNTAP  (see  [12], 
Theorem  3.16  and  Corollary  3.7),  we  prove 

II^-1(2/)IIa?(x„)  < ci||F-1(2/)||(Xo,x1)x^.g.  (10) 

l|f-1(2/)ll(x0,xi)^„  <C2||F-1(3/)|U?(Xo)-  (11) 

By  obtaining  appropriate  upper  bounds  for  the  K-functionals  in  the  definition 
of  (7f0,  Xi)giq  and  (Fo,  FjJ^,,  0 < 9 < 1,  and  using  the  embeddings  X\  > X0, 
Fi  <— > Po,  it  can  be  shown  that,  for  Lipschitz  operators  F,  F-1, 

l|ir'‘1(y)ll(Xo,x1)7i_,,  < c3(II-F1_1(o)|[x1  + IlyllcYb.yiiyj-,,))  (i2) 

< c4(||F(0)||n  + ||x||(XoiXl)^  J (13) 

hold.  Combining  (10)  with  (12)  and  (13)  with  (11),  we  arrive  at  (8)  and  (9), 
respectively.  □ 
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Corollary  1.  Under  the  conditions  of  Theorem  2,  let  0 < s < A - a,  and 
assume  that  F( 0)  = Oy,,  F_1(0)  = Ox!-  Then, 

(Yo,Y i)^,,  - Al(X0)  - (Yo.Yi)^,,.  (14) 

In  particular,  if  a = (3  — 0,  then 

(y0lYi)il,  = ii;(Xo)  (15) 

( isomorphism  of  the  spaces,  equivalence  of  the  Hilbert  norms). 

Note  that  this  special  case  corresponds  to  sublinear  operators. 

Remark.  If  the  dependence  of  C in  (4)  and/or  D in  (5)  on  N is  weaker  than 
polynomial,  e.g.,  logarithmic,  then  the  left-hand  and/or  right-hand  embedding 
in  (14)  can  be  sharpened  by  setting  a = 0 and/or  (3  = 0 and  modifying  the 
index  q.  We  omit  the  details. 

Multiresolution  Galerkin-Petrov  methods  for  monotone  operators  (Ex- 
ample 1)  are  included  as  partial  cases  in  Theorem  2 and  Corollary  1.  For 
monotone  operators,  we  have  Y0  = Xq(Hq),  Y)  = X{(Ifi),  where  H0,  Hi 
are  Hilbert  spaces  with  Hq  ' H\  which  are  sufficiently  far  away  from  each 
other  so  that  Xo  ’ X\  and  Xq{Hq)  Xl{H\)  hold  simultaneously.  Here 
X\  fl  Y\  is  assumed  to  be  dense  in  Hq  and  Hi.  The  projectors  P/v  and  Qn 
in  Example  1 are  assumed  generated  by  multiresolution,  which  ensures  that 
Gn  ~ Gn-  1 C G2N- 

In  the  rest  of  this  section  and  in  the  next  section  we  shall  discuss  how 
to  reduce  the  rates  a and  (3  in  Theorem  2 and  Corollary  1 to  zero  in  the 
presence  of  near-degeneracy.  To  this  end,  we  shall  study  the  analogue  of 
the  phenomenon  of  near-degeneracy  with  multiresolution  methods  based  on 
biorthogonal  wavelets. 

One  equivalent  norm  in  the  inhomogeneous  potential  spaces  Hs  (cf.,  e.g., 
[14,4]  for  p = q = 2)  is  given  by 


oo  2n  — 1 

II/IIh.  X { £ |a0*l2  + £22>s  £ £ l/?l'il2}1/2,  (16) 

j= 0 k £7Ln  ;=i 

with  0 < s < r,  where  in  [8]  r is  the  Lipschitz  regularity  of  the  compactly 
supported  scaling  functions  ip  € HT , Cp  € Hr  and  wavelets  € Hr , G 
Hr  of  the  biorthonormal  wavelet  bases,  with  respect  to  which  / G IIs  can  be 
expanded  as  follows: 


°o  2"-l 

f(x)  ~ £ a0k<P0k{x)  + £ £ £ a.  e.  x,  (17) 

k£7Ln  j=0k£7Ln  1=1 

where  a0k  = ( f,<fiok)L2 , Pjl  = Each  hypercube  in  the  Calderon- 

Zygmund  decomposition  of  lRn  and  Stein’s  construction  of  Whitney-type 
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extension  operators  corresponds  to  2"  — 1 basis  functions  , l = 

1,  ...,2”  — 1,  in  each  of  the  two  biorthonormal  bases.  The  convergence  in 
(17)  is  in  the  norm  topology  of  Hs,  but  also  Lebesgue  a.e.  on  the  domain  ft 
of  the  functions,  ft  may  be  ]Rn,  hyperrectangle,  correspond  to  the  periodic 
case,  or  even  general  Lipschitz-graph  domain.  We  refer  to  the  currently  most 
advanced  work  on  this  topic  [3] , as  well  as  to  the  extensive  account  [4]  (for  the 
case  of  homogeneous  potential  spaces,  see  [7],  in  the  special  case  p = q = 2). 

Definition  5.  Let  ji  € IN.  A non-degenerate  wavelet-based  projector  (NWP) 
is  denoted  by  Pj1  and  defined  by 

ji-l  2"-l 

pjJix)  = Ya°k^ok^  + Y Y Y Pjl'l’jli*)’  * ^ ft,  (18) 

k j—0  k (=1 

cf.  (17).  A near-degenerate  wavelet-based  projector  (NDWP)  is  denoted  by  Ph 
and  defined  by 

_ J(jl,k)-l2n-l 

Pjj(x)  = '52°‘ok<pok(x)+J2  Y Y x e (lQ) 

k k j= 0 1=1 

Vk  J(ji,k)  > J(ji  — 1,  ft),  J(ji,k)  >ji;  3kh  : J(ji,kh)  = ji-  (20) 

In  other  words,  for  a NWP  J(ji,k)  = ji  = const,  uniformly  in  k.  Thus, 
NDWP’s  are  a specific  partial  case  of  lacunary  wavelet-based  projectors  (see 
the  concluding  remarks  in  [4],  subsection  6.2),  lacunarity  being  with  respect 
to  the  NWP  corresponding  to  Ji  :=  max*,  J(ji,k). 

Example  2.  One  example  when  near-degenerate  FEM  or  lacunary  wavelet- 
based  projectors  of  NDWP  type  are  needed  is  in  the  error  analysis  of  numerical 
solutions  in  the  immediate  neighbourhood  of  the  boundary  dft  (see,  e.g.,  [11], 
Fig.  3.14,  3.15,  6.14,  6.15,  8.12).  Then  it  is  desirable  to  ensure  that  the  local 
approximation  rates  near  and  on  dfl  do  not  deteriorate  compared  to  the  local 
approximation  rates  in  the  interior  of  ft.  Indeed,  assume  that  5ft  is  regular 
enough  (Lipschitz  or  smoother).  Then,  by  the  trace  theorem  (see,  e.g.,  [1,9]), 
if  / 6 Hs(Ll),  ft  C 1R",  then  the  restriction  of  / on  5ft  is  less  regular,  namely, 
/|an  ^ /Ts~1/2(5ft)  holds.  Then,  the  local  approximation  rate  achieved  via 
NWP,  given  in  (18),  is  0( 2-A«)  jn  the  interior  of  ft  and  only  0(2-J1(s-1/2)) 
near  5ft.  To  achieve  the  desired  uniform  distribution  of  the  error  in  the  interior 
and  near  the  boundary  when  / is  smooth  enough  ( s > 1/2),  NDWP  given  in 
(19,20)  should  be  employed,  with  J(ji,k ) « j\  for  k corresponding  to  the 
interior  of  ft,  and  with  J(j\,k)  x Ci ji  + C2  otherwise,  where 

C2  >0,  <7i  = 1+2(s-1/2)  >L  (21) 

In  the  context  of  Theorem  2 and  Corollary  1,  if  Xj  = HSj , j = 0, 1,  with 
Sq  < si  so  that  X\  Xq  is  fulfilled,  then  it  can  be  verified  that  for  NWP 
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both  the  direct  inequality  ||/  - Pjlf\\x0  < Ci2  -^UfUxt  and  the  inverse 
inequality  ||f>j1/||A’i  < C,22-JlA||-f,j1/||x0  hold,  with  A = si  - so  and  with 
Ci,  C2  independent  of  j\.  Hence,  in  Definition  4 a = (3  = 0 is  attained.  On 
the  contrary,  for  NDWP  satisfying  (21)  the  constants  C\  and  C2  depend  on 
ji  and  a > 0,  /3  > 0 holds. 

In  the  case  of  NDWP,  is  it  possible  to  somehow  reduce  a and  /3  to  zero, 
thereby  achieving  isomorphism  in  (14)?  It  turns  out  that  the  answer  is  posi- 
tive, and  below  we  shall  propose  a general  method  how  to  achieve  this. 

Our  approach  will  be  consider  more  general  spaces  Xo,  X\  than  IIs , so 
that,  for  the  new  Xo  and  Xi,  a = (3  = 0 holds.  Consider  the  Hilbert  space 
Hs'w  with  norm 

00  2"  — 1 

ll/ll*-  *(£h>*l2  + ££22u'<^5  £ I/?Si2)1/2- 

k j= 0 k 1=1 

The  spaces  from  this  scale  still  admit  atomic  decomposition  via  the  same 
Riesz  bases  of  biorthonormal  wavelets  as  Hs . The  weight  w(j,k)  is  positive, 
monotonously  increasing  function  in  j for  each  fixed  k,  and  depends  on  the 
choice  of  J{ji,k)  in  (20).  The  definition  of  w(j,k)  is 

w(J(ji,k),k)=ji,  (22) 

w(j,k)=ji-l,  j = J(ji  - l,k),  J(ji  - l,fc)  + 1,  ...,  J(ji,k)  - 1,  (23) 
Vji  € IN  Vfc  € ZZn. 

Now,  take  Xj  = HSj’w,  j = 0,1,  with  sq  < s 1.  It  can  be  seen  that 
Xi  t— ► Xo  holds,  and  we  can  consider  this  pair  of  spaces  in  the  context  of 
Theorem  2 and  Corollary  1. 

Corollary  2.  Under  the  conditions  of  Corollary  1,  assume  that  Xj  = HSj'w , 
j = 0, 1,  where  w = w(j,  k ) is  the  left  inverse  ( see  (22,23))  of  J(j,  k)  as  defined 
in  (20).  Assume  also  that  N = 2J1  and  Gn  = PnX0,  where  the  NDWP  Pn 
is  defined  in  (19),  with  the  same  J(j,k)  in  (20).  Let  s : 0 < s < A = sj  — So- 
Then  (15)  holds. 

Proof:  (Outline.)  It  can  be  verified  that  the  bounds 

OO 

II/--p*/IIh-o.-<(£  £ £1/351 12)1/2>  (24) 

j=j  1 k:w(j,k)>j 1 I 

114/11*—  < (£|ao*|2  + £ £ 22^-fc)Sl£|/3j']|2)1/2,  (25) 

k j= 0 k:w{jyk)<j\  l 

hold.  (Recall  that  j\  = min^  J(j\,k),  J\  = max*,  J(j\,k).)  After  some  com- 
putations, (24)  and  (25)  imply 


||/-4/IU„<Ci2-^||/||x1,  V/eXi, 


(26) 
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II4/IU,  < C22^\\Phf\\Xa,  if  e V0>  (27) 

with  A = Si  — so,  and  the  constants  C\  and  C2  in  (26,27)  do  not  depend  on 
j 1,  i.e.,  for  this  choice  of  the  spaces  Xq,  X\  in  Definition  4 » = /3  = 0 holds. 
The  result  now  follows  from  Corollary  1.  □ 

Thus,  we  have  solved  the  problem  of  characterizing  the  best  approxima- 
tion spaces  induced  by  NDWP  defined  in  (19)  and  (20).  In  this  approach  we 
remained  entirely  within  the  classical  BNTAP.  There  is  also  another  approach 
which  goes  beyond  the  general  BNTAP,  by  abandoning  the  use  of  the  real 
interpolation  functor.  This  approach  leads  to  atomic  decomposition  of  Wiener 
amalgam  spaces  and  will  be  considered  elsewhere. 
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